Lorentz Symmetry Violation and Galactic Magnetism 
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We analyze the generation of primordial magnetic fields during de Sitter inflation in a Lorentz- 
violating theory of Electrodynamics containing a Chern-Simons term which couples the photon to 
an external four-vector. We find that, for appropriate magnitude of the four-vector, the generated 
field is maximally helical and, through an inverse cascade caused by turbulence of primeval plasma, 
reaches at the time of protogalactic collapse an intensity and correlation length such as to directly 
explain galactic magnetism. 



Up to today, no compelling evidence for violation of 
Lorentz symmetry exists (for a recent review see 
The main motivation to investigate possible consequences 
of Lorentz violation is that in unified theories at the 
Planck scale, such as quantum gravity, it is allowed or 
even unavoidable. Hence, any direct or indirect evidence 
of Lorentz violation would give information about very- 
high-energy physics or, in a cosmological context, about 
physics of the very early universe. 

In this paper, we indeed claim that the existence of 
galactic magnetic fields could be an indirect evidence 
of Lorentz symmetry violation. Magnetic fields are ob- 
served in any types of galaxies (for reviews see jH). They 
are correlated on scales of order of galactic dimensions 
(~ lOkpc) and have an intensity of about 10~^G. 

In the existing literature, it is assumed that a pri- 
mordial magnetic field (coherent within domains of 
physical size Cphys) remains frozen into the plasma dur- 
ing the collapse of a protogalaxy owing to the high 
conductivity of the protogalactic gas 0. Accordingly, 
since the physical size iphys — lOkpc and mass density 
Ppg ~ 10~^'^gcm~^ of a protogalaxy remain approxi- 
mately constant during the collapse, is amplified by a 
factor [ppsl pit)Y^^ while ^phys is reduced by [ppg/ p(t)Y^^ 
with respect to a field which is not trapped in a proto- 
galaxy. Here, p{t) oc \/ is the (time-dependent) total 
energy density of the universe. Taking into account that 
"^^gcm"-^ 



p{h) - 10 
tual age of the universe 



where ~ 13.7 Gyr is the ac- 
, one find the "canonical" result 
that the overall amplification of B^ and diluition of i^phys 
are respectively 10'* and 10^ [2]. 

The above analysis, however, is not completely cor- 
rect. In fact, being the collapse of the protogalactic gas 
a turbulent phenomenon [6^ , we expect that the magnetic 
field will be processed by magnetohydrodynamic effects. 
Defining as usual the kinetic and magnetic Reynolds 
numbers at the length scale ^phys as Re = w^phys/^^ 
and Rcb = w^physCc where v is the typical velocity of 
the fluid motion, v the kinematic viscosity, and dc the 
conductivity, magnetohydrodynamic turbulence occurs 
when Re and Res are greater than unity Apply- 
ing the results of 0] we have that, in a typical proto- 
galaxy, V is given by the "virial thermal speed" , Wth — 
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the kinematic viscosity is dominated by ions, 
1/ ~ z/; ~ lO-'^pc, while the conductivity is the inverse of 
the Spitzer resistivity, rjs — lO^^^pc. Consequently, we 
have Re ~ ^phys/O.lpc and Res ~ lO^^Re, which shows 
that magnetic domains of physical size greater than O.lpc 
evolve turbulently during protogalactic collapse. 

In the framework of freely-decaying magnetohydro- 
dynamic turbulence, it is a well-known result that the 
evolution of a maximally helical field is characterized 
by the quasi-conservation of magnetic helicity which, in 
turn, triggers and sustains a transfer of magnetic en- 
ergy from small to large scales, a mechanism known as 
inverse cascade (for a definition of magnetic helic- 
ity and maximally helical fields, see below). As a re- 
sult, the intensity of the magnetic field inside a domain 
is not longer constant but decays as B^ cx Cph/s^ 
while the physical correlation length grows as ^phys(i) — 
fphys(ipg)[i/7"oddy(ipg)]^^^ 9] for t > Tcddy(ipg), whcrc 
tpg is the time when protogalactic collapse begins (cor- 
responding to a redshift of order Zpg ~ 10^ [2., 10]) and 
Toddy (i) — fphys/w is the so-callcd eddy turnover time. 
Using this latter evolution law, we find that after a time 
~ [0.1pc/^phys(ipg)]*^^io the correlation length reaches 
the maximum possible value corresponding to the pro- 
togalaxy dimension iphys- This means that for t > t^, 
the magnetic field remains effectively frozen in the pro- 
togalaxy, B^{t) = const and ^phys(0 = const. The fore- 
going analysis, if on the one hand is oversimplified, on 
the other hand indicates that if the magnetic field at 
the time of protogalactic collapse is i) turbulent [i.e. 
^phys(ipg) ^ O.lpc], a) maximally helical, and in) has an 
intensity equal to ^^(ipg) - 10"^[10kpc/^phys(ipg)]*/^G 
[so that B^{t^,) — lO^^G], it explains galactic magnetism. 

The possibility that Lorentz violation could give rise to 
large-scale magnetic fields was investigated for the first 
time by Bertolami and Mota [ll| in the context of string 
theories. Developin g th e original idea by Kostelecky, 
Potting and Samuel [ij], they showed that the appear- 
ance of an effective photon mass owing to spontaneous 
breaking of Lorentz invariance allows the generation of 
astrophysically interesting magnetic fields within infla- 
tionary scenarios. To our knowledge, the connection be- 
tween Lorentz symmetry violation and cosmic magnetic 
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fields has been studied only in three other papers, within 
the framework of either noncommutative spacetimes [isj 
or quantum theory with noncommutative fields [l4 1 . 

In this paper, we consider a modification of standard 
electromagnetism by adding to the Maxwell Lagrangian 
a Chern-Simons term. Although this is not the only ad- 
missible Lorentz- violating term [15], it is the archetypal 
one. First studied by Carroll, Field, and Jackiw \2£^ in 
1990, its cosmological and astrophysical implications are 
sill investigated nowadays (see, e.g., Q)- In curved space- 
times, the Maxwell-Chern-Simons action reads llSl: 



Sues = Jd^xel^-l F^^F^'' - \ v^e^A.F^^ 



(1) 



where F^,^ — d^A^, — d^A^ is the electromagnetic field 
strength tensor, F^''' = {l/2e) e^'^'P" Fp„, e the determi- 
nant of the vierbein e^, and e^^'^f" the Levi-Civita tensor 
density. The external four- vector p'' pick out a preferred 
direction in spacetime and then breaks (particle) Lorentz 
invariance 3 12 1 Si- Moreover, the Chern-Simons term 
violates CPT symmetry and, in general, conformal in- 
variance (which is an essential requirement for the gen- 
eration of seed magnetic fields in spatially- flat models of 
the universe [l9jV However, if one wants to retain gauge 
invariance, a supplementary condition on has to be 
imposed [la|: D^^p^ - D^p^ = d^py - d^p^j, = 0, where 
is the spacetime covariant derivative. In particular, 
the above condition is satisfied if is a constant, or if 
= and p'^ is an arbitrary function of time (Latin 
indexes run from 1 to 3, while Greek ones from to 3). 

The equations of motion follow from action ([1]), 
F>^F'^'^ _ p^pfj^:^^ while the Bianchi identities are 
jj^ptj-v — 0. We assume that during inflation the 
universe is described by a de Sitter metric ds^ = 
a'^{drf' — dx^), where a{ri) is the expansion parameter, 
r] = —l/{aH) the conformal time, and H the Hubble pa- 
rameter. Let us introduce in the usual way the electric 
and magnetic fields as F^i = —a^Ei and Fij = eijkO^Bk- 
In vector notation, the equations of motion read 



and inflation-produced magnetic fields will appear ho- 
mogeneous and isotropic. This means that the two-point 
correlation tensor of the magnetic field, (i?i(x)i3j(y)), 
where (...) denotes ensemble average, is a function of 
|x — y| only and moreover transforms as an S'0(3) ten- 
sor. In terms of the Fourier amplitudes of the magnetic 
field, Bk = /cJ^x e"*'' '' B(x), where k is the comoving 
wavenumber, these conditions translate into 21|: 



((Bk).(BkO,) ='5(k + k') 



p(S)7lf^ _ (A)7lfHfe 
ij ^3 ij ^2 



, (3) 



where P/f^ = {2'Kf{6,, - hk 



,(A) 



1] ^ V""/ i^ii^jji ij — i{2Tr)^Sijiki, Sij 



is the Kronecker delta, ki = ki/k, k — |k|, and Sijk is 
the totally antisymmetric tensor in tree dimensions. The 
functions Vk and Hk, the so-called magnetic field and 
magnetic helicity power spectra, denote the symmetric 
and antisymmetric parts of the correlator. 

Let us define the average magnetic field on a comoving 
scale A and the magnetic helicity density in a volume V = 
respectively a.s ^ Bf = ( | /d^y H/^(x - y) B(y)|2) 
andHx ^ { \ Jd^yJd^zWxix-y)Wx{:K-z) A{y)-B{z)\), 
where Wa(x) = {2TrX^)~^^^e~^^^ ^^^^ ' is a Gaussian win- 
dow function. Then, taking into account Eq. we get 

Bl= j^^Wl{k)Vk, Hx = J^^fwrn^k, (4) 

with Wx{k) — e^^ the Fourier transform of Wa(x). 
It is useful to introduce the so-called orthonormal he- 
licity basis {ej,e^,e^}, where = — i(ej^ ± «e^)/-\/2, 
= k, and {6^,6^, e^} form a right-handed orthonor- 
mal basis. The Fourier transform of the magnetic field 
can be then decomposed as Bk — B^e^ + B^e^, B^ 
being the positive and negative helicity components of 

Bk, so the magnetic and helicity power spectra be- 
I23I1 15, — /-^,3/o^2^ 



come 



Vk = {k'/2n') 



\B^\^) and TLk 



5,,(a2E) - V X a^B 



ma^B 



p X a E 



(2) 



(fc^/27r^) (|i?k 1^ ^ l-^k P)- ^ magnetic field is maximally 
helical at the time 77' if either B^[rj') or B^{ri') is zero. 

Taking the curl of Eq. ([2]), using the Bianchi identities, 
and going into Fourier space, we get 



and V • E = — p • B, while the Bianchi identities become 
dr^{a^'B) + Vy.a^'E = and V-B = 0, where = (to,p). 

From Eq. ([2]), we readily get that the effect of p 
is just to cause a precession of the electric field vec- 
tor. Therefore, if the external four-vector takes the form 
pP — (0,p), the Chern-Simons term does not break con- 
formal invariance of electromagnetism so that, as shown 
long time ago by Turner and Widrow , no astrophys- 
ically interesting magnetic fields can be produced during 
inflation. For this reason, we will restrict our analysis to 
the case p^ = (m, 0), with m a positive mass parameter 
which, eventually, depends on the time. In this case, the 
external four-vector does not introduce any preferred di- 
rection in space (but just a preferred direction in time) 



dlia^B^) + k(k ± m) a^B^ = 0. 



(5) 



Before proceeding farther, we note that Eq. ([5]) is valid 
also during radiation and matter dominated eras. Posit- 
ing in it the phase-exponential ansatz, a^B^ cx e*"'', 
we recover the well-known result that two polarization 
states of photons propagate with different phase veloci- 
ties, causing the polarization plane to rotate by an angle 
A0 — ~{l/2)Jdr]m{r]) (to the first order in m/k) [l6| . 

The study of this effect on cosmological scales, known 
as cosmological birefringence, gives the most stringent 
bounds on m 0, [2J]- In fact, Ac/) has been constrained 
by observation of radio galaxies and quasars, A0(z3C9) = 
2° ±3° (68% C.L.) [H (where 3C9 at Z3C9 = 2.012 is the 
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best data sets available at a single redshift), and by the 
analysis of Cosmic Microwave Background (CMB) polar- 
ization fluctuations In the latter case, A;/) has been 
constrained between the reionization epoch (zroion ^ 10) 
and today, -22.2° < A(/i(zreion) < 7.2° (95% C.L.), and 
between the decoupling epoch (zdoc — 1090) and today, 
-5.5° < A0(zdoc) < 3.1° (95% C.L.). The above limits 
give, with m constant in time, the bound m < 10~^*eV. 
However, this bound refers to relatively recent epochs 
and it is not excluded that in the early universe the 
value of m could have been much greater. In order to 
avoid inessential complications, we assume that Lorentz 
symmetry is broken at some energy scale before or dur- 
ing inflation, and that m is constant during inflation and 
decreases afterward to fulfil the above constraint. 
The solution of Eq. ([5]) with m constant is 



a^B^ = ±kexp[ikrj^/T±~m/k]/V2k. 



(6) 



Here, the normalization is such that for j] — oo (where 
m is supposed to be zero) we have the usual Maxwell 
vacuum spectrum, a^B^ = zLke^''^ /V2k. From Eq. ^ 
we easily deduce that only negative-helicity modes such 
that 1/ (mrf') <^ k <m are (exponentially) amplified and 
that in this case a^B^ ~ kexp{k\ri\y/m/k)/^/2k. Using 
Eq. (|3|), the amplitude of the (almost) maximally helical 
field at the end of inflation, rj = yyrh, reads 



a,\BA(?7rh)=^eV27rA2V^, 



(7) 



where aj-h — a(?7rh), 7 — ("^'7rh/'^)"'^^^' assumed 
1/"^ ^ ''TT-Vrh 126]- Looking at Eq. ([7]), we see that 
the magnetic field is peaked at small scales since the ex- 
ponential factor is inversely proportional to ^/X- There- 
fore, at the end of inflation, magnetic power is highly 
concentrated in scales just a little greater than A ~ 1 /m 
and negligible at larger scales. In this case, the condition 
A ^ TO'yrh (or equivalently 7 1) reads m|77rh| ^ 1- 
This means that the mode corresponding to A ~ 1/m is 
well inside the horizon since its physical wavelength is 
much smaller than the Hubble radius H~^. 

After inflation, the universe enters in the so-called re- 
heating phase, during which the energy of the inflaton is 
converted into ordinary matter. In this paper, we restrict 
our analysis to the case of instant reheating, that is we 
assume that after inflation the universe enters the radia- 
tion dominated era. In this era, as well as in the subse- 
quent matter era, the effects of the conducting primordial 
plasma are important when studying the evolution of a 
magnetic field. They are taken into account by adding to 
the Maxwell action the source term Jd^xejM.^, where 
— (0, (TcE) is the external current [2^ (we are as- 
suming that the Chern-Simons term is negligible after 
inflation). Plasma effects introduce, in the left-hand-side 
of Eq. (O, the term —aacdr^{a'^B^). 

For modes well inside the horizon, fc|?7| ^ 1, Eq. ^ 
reduces to k'^a'^B^ ~ —aacdrf{a'^B^) so that a?B^{ri) = 



"rh^k ('^rh) exp[— J'd77(fc^/a(Tc)]- Introducing the comov- 
ing dissipation length, = j dri(l/ aac) — l/ia^Hoc), 
where we used the relation 77 ~ 1/ {aH) valid in all eras, 
we get a^B^ijf) = a\B^{T^,h) exp(-^^i3,fc2). This means 
that modes with comoving wavelength greater than the 
comoving dissipation length are frozen into the plasma, 
while modes for which A < ^diss are washed out. In par- 
ticular, only modes with comoving wavelength greater 
than the actual dissipation length, ^diss(7o) — lAU [1], 
can survive until today (Tq ~ 2 x 10~^eV is the actual 
temperature But this is not our case since wave- 

lengths greater than lAU would correspond to modes 
outside the horizon at the end of inflation. Thus, we find 
the uninteresting result that the inflation-produced field 
is dissipated soon after the beginning of radiation era. 

This situation drastically changes if one takes into ac- 
count that the primordial plasma could have been turbu- 
lent [3] • Indeed, turbulence naturally arises during (first- 
order) cosmological phase transitions and those could 
have affected the universe during its evolution. In the 
following we simply assume that the early universe is tur- 
bulent without referring to any specific, generating model 
of turbulence (if this is not the case, our mechanism can- 
not account for the presence of large-scale magnetic fields 
in galaxies, as discussed above). 

At the beginning of radiation era, T — Trh, and at the 
comoving scale 1/m the kinetic and magnetic Reynolds 
numbers are Re ~ avT^h/m and Res ~ vTr\i/{am), 
where we used ~ 1/ (aT) , CTc — T/a Q , and a is 
fine structure constant. For an incompressible, relativis- 
tic turbulent plasma an upper limit to v is the sound 
speed Us = l/\/3. Taking, for example, Trh — lO^GeV, 
m ~ lO^^^eV, and v ~ Vs, we get Re ~ 10^^ and 
Res ~ 10'^^. The temperature at the beginning of radia- 
tion era, the so-called reheat temperature, is constrained 
in the range IGeV < Trh ^ lO^GeV from analyses of Big 
Bang Nucleosynthesis (BBN) and CMB radiation 27 1. 
For instant reheating, it results Trh — {Hmpiy^^ 20|, 
where H is the Hubble parameter during inflation and 
TOpi ~ lO^^GeV the Planck mass. 

In a spatially-flat Friedman universe, the decay of a 
turbulent, maximally helical field proceeds through in- 
verse cascade and the evolution laws are equal to those in 
a static Minkowski universe provided that time, intensity 
and physical correlation length of the magnetic field are 
replaced by the conformal time, the "comoving intensity" 
a^B(^ and comoving length, respectively Therefore, 
we have a'^B^ oc ^"^/^ and ^(77) ~ ^(r?rh)[??/Teddy(r?rh)]^''^, 
where Teddy(?7ih) = 'C('7rh)/w(?7rh)- The inverse cascade 
terminates when the primordial plasma ceases to be tur- 
bulent. This happens when the universe has sufliciently 
cooled down to render kinematic viscosity very large, 
thereby causing the kinetic Reynolds number to drop be- 
low unity at the scale of interest. From this point on, 
the evolution of the magnetic fleld is adiabatic, in the 
sense that a^B^ = const and ^ = const. Taking into 
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account that for T <C IMeV the kinematic viscosity is 
V ~ lO^^MeV^ T'^ M ai^d that v scales in time as the co- 
moving magnetic field, v{ri) ~ w (/?rh) [Teddy ('yrh)/'^]^/^ 0, 
we obtain the temperature at which the plasma becomes 
non-turbulent [starting at the time ?7rh with w(??rh) — '^s 
and C(77rh) ~ 1/m]: T^t ~ i2m%^eV or Tnt 28mf/"eV 
according to mi3 ^ 6 x 10^^ or mi3 > 6 x 10^^, where 
mi3 = m/10~^'^eV. The two cases correspond, respec- 
tively, to have T„t > Teq or Tnt < Teq, where Teq ~ 3eV 
is the temperature at the matter-radiation equality [3| ■ 
Evolving (along the lines discussed above) the mag- 
netic field ([7]) evaluated at the scale 1/m, from the end 
of inflation till the time of protogalactic collapse, we get 

-,1/2 



B^(tpg) ~ 3 X 10 

fphys(ipg) - 



lOkpc 

^phys (^pg 
-1/2 



G, 



0.1 TOj^o pc if mi3 > 6 X 10 



0.3 TO 



-7/17 
13 



pc if mi3 < 6 X 10 



(8) 
(9) 



where 7 ~ m|r/,h| ~ ^2mi3T~^ and Tg = T,h/10^GeV. 
In order to explain the existence of galactic magnetic 
fields, we impose Bj(ipg) - 10~^[10kpc/^phys(ipg)]^/^G. 
From Eq. ([5]) , it is easy to see that the above requirement 
is satisfied for mi3 ~ Tg, so we find that the magnetic 
field at the time of protogalactic collapse has a physical 
correlation length ranging from O.lpc to 0.6kpc and then 
could be exactly the field detected in galaxies [llj. 

Finally, if we assume as in Ref. ll| that to scales as a 
power of temperature, m{T) oc T', and take I = 1 (i.e. to 
simply redshifts after inflation), we find that its present- 
day value is too ~ 2 x 10~'^'*eV, and that the rotation an- 
gle of polarized light emitted from a source at a redshift Zi 
and detected today is A0(zi) ~ (3/2)ioTOo[l — ^/TA~Zi]. 
For the three cases discussed above, we get A(/)(z3C9) — 
-8°, A(/.(2,eion) ^ -26°, and A<^(2:dec) ^ -362° (equiv- 
alent to —2°) which shows, roughly speaking, that m is 
consistent with observational constraints for ^ ^ 1. 

It is worth stressing that, although other mechanisms 
for generating cosmic magnetic fields exist in the litera- 
ture [Ij, the field produced in our mechanism has three 
features which are not simultaneously present in any of 
them: it directly explains galactic magnetism (without 
resorting to further and still-debated amplification mech- 
anisms such as galactic dynamos 0, [l^ , it is maximally 
helical, and its power is concentrated on a single scale, 
the magnetic correlation length. These characteristics 
make it peculiar and hopefully detectable in future ob- 
servations of CMB temperature and polarization fluctu- 
ations (30| . 

In conclusion, we have studied the generation of a 
cosmic magnetic field in Maxwell-Chern-Simons theory 
of Electrodynamics containing a Lorentz-violating cou- 
pling between the photon and an external four-vector. 
The breaking of conformal invariance in this theory al- 
low electromagnetic vacuum fiuctuations to be ampli- 
fied during inflation. The resulting field at the end of 



(de Sitter) inflation is maximally helical and is subse- 
quently processed by magnetohydrodynamic effects tak- 
ing place in the primeval turbulent plasma. In particular, 
quasi-conservation of magnetic helicity is responsible for 
an inverse cascade of magnetic spectrum which causes a 
growth of the correlation length together with a slow dis- 
sipation of the magnetic field intensity. As a result, for 
appropriate magnitude of the four-vector, the produced 
field at the time of protogalactic collapse has an intensity 
and correlation length of the right magnitude in order to 
explain galactic magnetism. 

We thank V. A. Kostelecky for valuable comments. 
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